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Comment on "The third Zemach moment of the proton", by Cloet and Miller. 
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Cloet and Miller, in arXiv: 1008.4345, state that "existing data rule out a value of the third Zemach 
moment large enough to explain the current puzzle with the proton charge radius determined from 
the Lamb shift in muonic Hydrogen. This is in contrast with the recent claim of De Rujula" . To be 
more precise: it is not. It is, however, contrary to what they claim that I claim. Cloet and Miller 
have simply misinterpreted my claims pp. 

PACS numbers: 31.30.jr, 12.20.-m, 32.30.-r, 21.10.Ft 



I. MY MAIN CLAIM 



II. MY ALLEGED CLAIM 



Let "model-independent" mean "independent of a par- 
ticular parametrization of the charge distribution of the 
proton, p p (r)". The most precise model-independent 
measurements of the corresponding mean square radius, 
(r 2 ), are mainly based on the theory [2] and observations 
[3J of Hydrogen. The result, compiled in CODATA @], is 



(r 2 ) (COD ATA) = (0.8768 ± 0.0069 hnf 



(1) 



Consider next the 2P^ 2 —> 2Sf^ 1 Lamb shift in a pp 
atom. In meV units for energy and fermi units for the 
radii, the predicted value [5] is of the form 

L th [(r 2 p ),(4)(2)] = 

209.9779(49) - 5.2262 {r 2 p ) + 0.00913 (r|) (2 ), (2) 

which is also model- independent. The first two coeffi- 
cients are best estimates of many contributions [B] while 
the third stems from the n = 2 value of an addend [7] 



AE 3 (n) = 3 m r Si (r p )(2), 
proportional to the third Zemach moment 



< r p>(2) 



d 3 r 1 d 3 r 2 p(ri)p(r 2 )\ri - r 2 | 



(3) 



(4) 



The quoted Lamb shift has been measured [6 to be 



-'cxp 



206.2949 ± 0.0032 meV. 



(5) 



Assume that the the ory a nd experiments quoted so far 
are correct. Use Eqs. ( 2|5 1 to write L th = L cxp . Intro- 
duce into the expression for L th the observed value of 
(r 2 ) given in Eq. (flj). Solve for (r p )( 2 )- The result is: 



pi 



WD 



(2)J 



I 1/3 



3.32 ± 0.22 fm 



(6) 



with the error dominated by the CODATA uncertainty. 

My main claim is to have done the algebra leading to 
Eq. J6I), which is model- independent. 



Cloet and Miller state "In the simple monopole model 
of De Rujula, which claims to account for the proton 
radius puzzle..." . 

I have not claimed that a monopole model (which, in- 
cidentally, is not "mine") accounts for the proton radius 
puzzle. On the contrary, I explicitly constructed a "toy 
model" with two added monolope form factors to illus- 
trate how the non-vanishing contributions of no fewer 
than two monopoles were needed to reconcile the data. 

The comment that may have led Cloet and Miller to 
misinterpret me is the last paragraph in arXiv: 1008.3861 

"The result Eq. ^ is ^(^-independent and to be 
treated with due respect. Right after offering excuses, 
I shall break this rule. The third Zemach moment is very 
sensitive to the long-distance part of p(r)... A correct 
question to ask is the scale, m, to which this tail corre- 
sponds. Not to be confused with the toy model fit to both 
(r 2 ) and (r 3 )^), the simplest naive -but consistent- an- 
swer is provided by Eq. ^ in the single-pole approxima- 
tion i.e. tn — 261 MeV, tantalizingly close to the threshold 
of the proton form factor's cut at 2m I ± ~278 MeV". 

All I did in the previous paragraph is to use a monopole 
to extract a rough estimate of the slope of p p (r) at large 
r, and to express it as an inverse mass. I did not contend 
that a single monopole representation of Ge{— q 2 ) "ac- 
counts for the proton radius puzzle" at all relevant values 
of q 2 , nor that its Fourier transform accurately describes 
p p {r) at all relevant r. 

The penultimate paragraph above, from v2 of 
arXiv: 1008.3861 (posted one day after vl) may have been 
more obscure in vl. No doubt this is what mislead Cloet 
and Miller, regarding what my claims really were. 



III. CLOET AND MILLER'S CLAIMS 

I use the word "claims" here for literary consistency. 
But I mean "statements" . 



2 



Cloet and Miller study the proton-radius puzzle by use 
of three models -the dipole and the ones in [3[9]- that de- 
scribe Ge(q 2 ) as measured in ep scattering experiments. 
They find that the model for which G^(9 2 ) falls fastest 
at large q 2 gives the largest value of (rp) 2 2 \/(r 2 ) 3 , but is 
insufficient to resolve the "puzzle" . 

In the models in [8, 9 the numerical parameters can 
be combined with their single mass scale (m p ) to be 
reinterpreted as different mass scales. This is in agree- 
ment with a claim that I have indeed made: "Any simple 
one-parameter description of the proton's non-relativistic 
Sacks form factor, Ge(— q 2 ), in terms of only one mass 
parameter is inaccurate: the proton is not so simple" . 

The crucial problem was adroitly emphasized by Sick 
10J. It is very difficult to extract reliable information on 
p(r) (such as its moments) from its Fourier transform, 
Ge(q) 2 . The radius of convergence of the expansion in 
r from which one directly extracts (r 2 ) from the ep data 
is so small, that one must use numerical simulations and 
a continued-fraction expansion to skirt the uncertainties 
associated with data normalization at small q 2 (and their 



systematic errors) and to obtain a stable, numerically- 
meaningful result not contaminated, for instance, by the 
term in (rp. Clearly, if extracting (r 2 ) is delicate, the 
more so it is to infer (r|)(2)- 

The authors of [9] are aware of the above-mentioned 
practicalities. But they do not need to face them a la 
Sick, since they analyze data at much larger q 2 than the 
range with which we are concerned. Cloet and Miller sim- 
ply extrapolate the high-q 2 models in [H [9] to the |q| val- 
ues of O (a mp) of interest here. In this way, not surpris- 
ingly, they obtain a model- and extrapolation-dependent 
value of (rp) 2 2 \ in disagreement with Eq. (K3|. 
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